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Abstract 

In spite of their original discrepancy, both dark energy and modified theory of gravity can be 
parameterized by the effective equation of state (EOS) oj for the expansion history of the Universe. 
A useful model independent approach to the EOS of them can be given by so-called Chevallier- 
Polarski-Linder (CPL) parametrization where two parameters of it (ujq and ui a ) can be constrained 
by the geometrical observations which suffer from degeneracies between models. The linear growth 
of large scale structure is usually used to remove these degeneracies. This growth can be described 
by the growth index parameter 7 and it can be parameterized by 70 + 7a (1 — a) m general. We 
use the scalar-tensor theories of gravity (STG) and show that the discernment between models is 
possible only when j a is not negligible. We show that the linear density perturbation of the matter 
component as a function of redshift severely constrains the viable subclasses of STG in terms of u 
and 7. From this method, we can rule out or prove the viable STG in future observations. When 
we use Z(<p) = 1, F shows the convex shape of evolution in a viable STG model. The viable STG 
models with Z(4>) = 1 are not distinguishable from dark energy models when we strongly limit 
the solar system constraint. 



1 Introduction 



The late-time acceleration of the cosmic expansion invokes the mysterious dark energy or the 
modification of the gravity theory beyond Einstein's general relativity. Although there exist 
a number of dark energy and modified gravity models, both can be described by the effective 
equation of state (EOS) uj when one considers the expansion history of the Universe. Chevallier- 
Polarski-Linder (CPL) parametrization uj = uq + co a (l — a) is one of the suitable candidates to 
describe the general uj [1, 2]. 

However, models of different physical origins with the same background expansion history can 
not be separated with uj. Thus, the growth of large scale structure is used as a complementary 
probe to segregate models [3, 4, 5, 6, 7, 8, 9, 10]. 

The growth index parameter (GIP) 7 defined by dh ^\™^ = ^ m («) 7 is often used to discrim- 
inate different models using the linear matter perturbations [11, 12, 13]. Although 7 is time 
dependent by its definition even for the simple dark energy models [14, 15, 16], the constant value 
of it can be well matched with some dark energy and modified gravity models [17, 18, 19, 20, 21, 22]. 
However, it can be generalized as a time dependent 7(a) [23, 24, 25, 27, 28, 29]. 

Behavior of 7 in the so-called Dvali-Gabadadze-Porrati (DGP) braneworld model [30, 31, 32] 
has been widely investigated [10, 19, 21, 26, 28, 33]. 7 for a modification of gravitational action with 
a general function of the scalar curvature instead of the standard Einstein-Hilbert term, named 
f(R) gravity [34] has been also studied [25, 35, 36, 37]. In general, 7 has a scale dependence in 
this model. Also when one considers the scalar-tensor theories of gravity (STG) [38, 39, 40], one 
can obtain the specific 7 for the certain STG model [24, 41, 42]. 

There is a useful approximate solution for the growth factor within general relativity [17] 

Sm (a) = ae f a W/«')lW-i] f (L1) 

where f2 m (a) is the matter density contrast, 5 m = 5p m /p m is the linear density perturbation of 
the matter, and 7 = 0.55 + 0.05[1 + oj(z = 1)] for uj > —1. However, the accuracy of this solution 
is misinformed. As either ojq or uj a increases, the accuracy of the approximate solution given in 
Eq. (1.1) is decreased. For example, the error of approximate solution is about 1.2% when one 
consider (ljq, u a ) = (—0.78, 0.32). If oj a increases to 0.4, then the error becomes 2.0% for the same 
value of wo- 

There are both theoretical and phenomenological motivations for STG. The former is related 
to the existence a ubiquitous fundamental scalar coupled to gravity in theories which unify gravity 
with other interactions [43, 44, 45, 46, 47]. Also the dynamical equivalence between f(R) theories 
and a particular class of STG has been shown in the case of metric formalism [48, 49, 50, 51, 52] 
as well as in the Palatini formalism [53, 54]. The later have several aspects. First, "the lithium 
problem" in the standard big bang nucleosynthesis (BBN) might be solved in STG due to the 
slower expansion than in general relativity before BBN, but faster during BBN [55, 56, 57]. The 
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weak lensing (WL) shear power spectrum in STG predicts the different one compared to in GR 
because they cause the different growth history of the matter [3, 58, 59]. Integrated Sachs- Wolfe 
(ISW) effect probes modified gravities on cosmological scales through the matter potential relation 
[60, 61, 62, 63]. The crossing phantom oj < — 1 also can be naturally obtained in STG [64, 65, 66]. 

There have been a number of reconstructions of specific STG models which is consistent with 
known observational constraints [67, 68, 69, 70, 71, 72, 73, 74, 75]. However, we need to reconstruct 
theory without any specific theory a priori. Thus, we use both background and growth history 
parameters (oj,j) to find the viable subclasses of STG. 

We briefly review the basic background evolution equations of STG model in the next section. 
In Sec. 3, we also review the linear perturbation equations of the model. We derive the recon- 
struction equations for model functions F{4>) and U(4>) as a function of scale factor a and check 
the viability of specific models in terms of parameters oj and 7 in Sec. 4. We conclude in Sec. 
5. We also show the accuracy of the approximate solution in Eq. (1.1) for the general values of 
(ojQ,oj a ) and find the initial values of 4>' and U in the appendix. 



2 Scalar-Tensor Theories of Gravity 

STG are described in the Jordan frame (JF) by the action [38, 39, 40] 

+ Sm^, V-m) , (2.1) 

where G* denotes the bare gravitational coupling constant which differs from the measured one, 
F(4>) and Z(4>) are dimensionless, and U{4>) is the potential of the scalar field 4>. F(4>) needs to 
be positive to ensure that the gravity is attractive. In the matter action S m (g^ Vl ip m ), the matter 
fields ip m is universally coupled to the metric g^ u and all experimental data including Hubble 
parameter H and redshift z will thus have their usual interpretation in this JF [76] . 

The evolution equations in the flat Friedmann-Lemaitre-Robertson- Walker (FLRW) metric are 
given by 

3FH 2 = 8ttG» (p m + p rad ) + X -Zft - 3HF + U , (2.2) 

2FH = -SttG* (p m + ^ p rad ) — Zcj) 2 — F + HF , (2.3) 

z(4> + 3Hj>) = ^F A R-^Z^ 2 -U A , (2.4) 
pi + 3H(l + 0Ji)pi = with (oj m = 0andw rac i = 1/3) , (2-5) 

where dots denote the differentiation with respect to (w.r.t) the (JF) cosmic time t, F^ = 
and R = 6(2H 2 + H) is the Ricci scalar. We limit our consideration after the matter dominated 
epoch and ignore the radiation component. The effective gravitational constant G c r between two 
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test masses measured in laboratory Cavendish- type experiments is given by [41, 77] 

G , ( 2Z^)F + 4F 2 A 
GcS -T[2Z^F + ^J- (2 - 6) 

This is different from the Newton's gravitational constant Gm = G*/F as the inverse factor of the 
curvature scalar R. In Eq. (2.6), G*/F comes from the exchange of a graviton between the two 
masses, whereas {G^/F)[F 2 ti /{2ZF + 3F^)] is due to the exchange of a scalar particle between 
them. When we use the familiar expression for Brans-Dicke (BD) representation, F = (f> and 
Z = ujbd/4>, Eq. (2.6) becomes Gf ff D = G^ - 1 (2lj B d + 4)/(2u BD + 3) = G^'^l + §/3 2 ) by 
using the conversion relation (3 = ' 4 ( 2tJ ^ D ^ gy between STG coupling constant and BD parameter 
[6]. We can recover G e s = G* in j3 — > (i.e. — > oo) limit. G e ff is experimentally bounded 
G cfr (z = 0) - Giv(z = 0) = 0.02% and G eff /G e// < 6 x lO" 12 ^- 1 [78, 79, 80]. One can always 
reduce Z((p) and F(4>) to one unknown function by a redefinition of the scalar field and thus we 
will consider Z(<fi) = 1 case. However, this parametrization Z(<p) = 1 can sometimes be singular 
[76]. 

There are several features we need to emphasize. First, the positive energy does not imply 
that (f> 2 > due to the mixture of tensor and scalar degrees of freedom in the JF. As the same 
reason, the second derivative of U does not give the precise value of its squared mass. Secondly, 
the evolution of the scalar field is determined by the effective potential U e g(4>) = U((f>) — \F{(f>)R 
as shown in Eq. (2.4). If we consider a light scalar field weakly coupled to matter (F ~ 0(1) 
and F $ ~ 0), then at early epoch U e s is dominated by F(<fi)R-term. Thus, the scalar field is 
dynamically driven to General Relativity corresponding value (4> = 0) and deviate from it at late 
epoch when U e g is approximated by U (</>) . 

In order to study the cosmological dynamics of the system, it is convenient to rewrite the 
evolution Eqs. (2.2)-(2.5) by using the new variable, n = In a and normalize the H 2 and H' by 
the present value of F 

3F H 2 = 8vrG*/9 m + ^H 2 <f)' 2 - 3H 2 F' + 3H 2 (F -F) + U, (2.7) 

- HH')F' + 2HH'(F - F) , (2.8) 

(2.9) 

p m + 3p m = 0, (2.10) 

where primes denote the differentiation w.r.t n, Fq means F{n = 0) and we only consider the 
matter component which is relevant to the late-time universe. We will use the above background 
equations (2.7) - (2.10) in followings. 
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3 Perturbations 



In this section, we review the matter density perturbations of STG in the longitudinal gauge w.r.t 
the JF [76, 81] 

ds 2 = -(1 + 2<5>)dt 2 + a 2 (l - 2^)dx 2 . (3.1) 



Also the energy momentum tensor components of the matter are given by 

T o = -(Pm + 5pm) , = -p m Vi 



(3.2) 



where p m (5p m ) is the energy density (contrast) of the pressureless matter and V{ is its velocity. 
The gauge invariant quantity is defined as <5 m = + 3Hv. The perturbed part of the energy 
momentum conservation in the Fourier space gives 
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k v 



, , ■ 3(* • Hr)'. 



a 2 H 



(3.3) 
(3.4) 



where k is a comoving wavenumber. On the other hand, the first order perturbed Einstein equa- 
tions give 



3F'$' + ^2A~ 2 F - Z<p' 2 + 3F'J $ 

$ _ $ 



8ttG*p 



jp -S m + I A" 2 -6- :;— \SF -;- — 



F /2' 

y 2 



#2 



F' 1 

+3— SF' + Z^'<5^ + 3Z<f>'8<l> + -<5Z0' 2 

_F 2 

87rG*p m ^: + Z</>'<50 + SF' - 5F , 

5F_ 

F ' 



(3.5) 
(3.6) 
(3.7) 



where A 



-i 2 H 2 



1 " " Also the perturbed part of the scalar field equation gives 

sr + {3 * + ^ + [ A - _ 3(2 + + _L(^) + ( 



ii 7 

A~ 2 ($ _ 2^) - 3(V + (4 + — )#' + 



From the equations (3.3) and (3.4), we obtain 



H 2 \ Z J,4> ' V Z 
^ + (3 *' + $')0 / -2|§|. 



C + ( 2 + ^)^m + A" 2 * = 3(¥ + + (6 + 3^) (vl/ + Hv)' . 
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(3.8) 



(3.9) 



We neglect time derivative terms w.r.t spatial derivative terms of corresponding perturbed vari- 
ables (i.e. subhorizon limit). This simplification holds at scales k > aH < lCT 3 /i/Mpc. In this 



4 



limit we can write the approximate equations for (3.5) and (3.8) 



A~ 2 $ ~ - — 
IF 



~ ($_ 2 ^)i^~ 9 $, 



(3.10) 
(3.11) 



where we use Eq. (3.7) in the second equality of Eq. (3.11). From the above two equations, we 
have the Poisson's equation 

k 2 

~ -47rG effPm <5 m , (3.12) 

cr 

where G c s is given in (2.6). Combining equations (3.9) and (3.12) gives 

< + (2 + §^ n -^^^0. (3.13) 

We can find the anisotropic parameter rj defined by rj = ($> — Vl/)/* [6] from the equations (3.7) 
and (3.11) 

F 2 , 



7] ~ — 



ZF + 3F 



■2 ■ 



(3.14) 



4 Constraints of F(<p) and U((f>) 

From the observational viewpoint, it is quite useful to use a parametrization of both the expansion 
history and the growth factor of the matter perturbation in terms of the EOS u)(n) and the GIP 
7(n). Chevallier-Polarski-Linder (CPL) parameterization uj = UQ+uj a {\ — e n ) uses its present value 
{ojq) and variation (uj a ) and it might be suitable for future observations [1, 2]. The growth factor 
can be given by 5 m (n) = e n g(n) where g(n) is given in Eq. (1.1) with 7 = 0.55 + 0.05(1 + uJdc{z = 
1)) [17, 23]. Although this functional form is accurate at certain level, there are several drawbacks 
in it. We explain the details about these problems in the appendix. 

Eq. (3.13) can be rewritten if we use the definition of the GIP 7 by using / = dln ^™ ( -"'- > = 
tt m (n)T [11, 12, 13] 

~2F + 4(F' /ct>') 2 ~ 



2F + 3(F> /ft) 2 



n m (n) 



«.<»)»- + (ym^w +7 |M + 2 + *y m[nV _ |£ 

=f f! m (n) 21 + (y lnf!m(») + 1~TT + 2 + ^WnP - ~!i„(n) =0, (4.1) 
V i'm(rcj ti J lb 

where we use the assumption |F(n)| S> \F^(n) 2 \ in the approximation. Thus, F(n) is given by 

F{n) 3 n m (n) 



F 2 P(n) 



, where 



P(n) = O m (n)^a n (n)^ + 7 / lnan(n)+ 7 §4 !! Y + 2 + ^ (4.2) 
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^m(n) 7 n m ( n y + i in n m (n) - 7 



3 + 2 



if 



+ 2 + 



where we use the relation O m (ri)' = —(3 + 2H' / H)Q m (n) in the second equality of -P(n). When 
we adopt CPL oj, then H 2 and H' /H are given by 
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= n m0 e' 3n + (1 - O m0 )e- 3 ( 1+W0+ ^) n e- 3 ^ (1 - en) 



1+co 



(1 - O m0 )e- 3 ^«+^) n e- 3 ^( 1 - en ) 



fi m0 + (1 - f] m0 )e- 3 (^»+^) n e- 3 ^( 1 - en ) 



(4.3) 
(4.4) 



Fq ~ 1 in order to satisfy G c ff(0) ~ Gat. After we obtain the functional form of F(n), we can also 
derive the that of the potential U(n) by using Eqs. (2.7) and (2.8) 



U(n) 



1H 2 F" 



F H 2 H Q \F 



+ 



5 + 
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IF 
H 



F 

y 



3J2 m (n) 



The effective equation of state uj is obtained from the same equations 

_ (f/ 2 + 2F" + 2(2 + H'/H)F' + 4{F - F )H'/H + 6{F - F ) - 2U/H 2 
^ ~ (j)' 2 - QF' - 6(F - F ) + 2U/H 2 

where the evolution equation of scalar field (j)' is given from Eq. (2.8) 



(4.5) 



(4.6) 



<}>' = \j~F" + (l-^)F'- 2^F - 3F a, 



(4.7) 



Thus, the EOS obtained from Eq. (4.6) should be same as that of CPL. As long as Fq ~ 1 is 
satisfied, the above constraint is well matched. We obtain the useful recursion relations for the 
differentiation of F from Eq. (4.2) 
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~F~ Q 



3 + 2- + 



P'- 


2 




rH'- 




rP'i 


T- 
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IH- 







\ F 



(4.8) 
(4.9) 



Now we probe the viabilities of some specific models which can mimic some well known dark 
energy and modified gravity models. We assume that 7 = 70 + 7 (1 — a) and $7 m o = 0.3 in the 
following analysis. One can always find the constant value of 7 (i.e. j a = 0) from Eq. (4.2) by 
using F(0)/F = 1. 



4.1 Comparison with Dark Energy models 

First, we investigate the models which mimic the cosmological constant (u)o,u) a ) = (—1.0,0) in 
the background evolution. We can further specify the models by considering the different growth 
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history (i.e. different value of 7). When 7 a = 0, 70 ^ 0.56 is found to get F(0) = 1. However, 
this gives (F y< p/\^F)\ = — 0.2 which violates the solar system test. The parameter set (70, 7a) = 
(0.56, -0.012) gives F(0) = 1.017 and (F >tj> /y/F)\o = -0.019 and passes the solar system test. But 
(f)' becomes imaginary at n < —0.6 (i.e. z > 0.8) in this case. This is quite similar to the result 
in Ref. [76] where F(ji)/Fq is defined as e~ 2n f(n) with f(n) is a general function of n. Thus, we 
may need to consider general Z(4>) in order to avoid this singular behavior of <fi' [82]. Even though 
the above case (0.56, —0.012) might be cured by introducing the nontrivial Z ((/>), it is difficult to 
be distinguished from the cosmological constant (A) model because the values of both oj and 7 are 
quite similar to those of A model. Thus, it is worth to check the case when 7 is quite different from 
that of A model. We try with (70,7a) — (0.6,0.126) which gives the quite interesting evolutions 
of physical quantities. Even though the potential U(n) is negative at present, it is fine because 
the second derivative of U(4>) does not give the precise value of its squared mass in JF. We show 
the evolutions of physical quantities in Fig. 1. As shown in the first column at the first row, oj 
obtained from Eq. (4.6) is exactly same as that of CPL with ojq = —1.0 and oj a = 0. Thus, if 
one investigate only the geometrical tests, then this STG theory can not be distinguished from A 

model. However, F(ti)/Fq can reach to 1.5 around n ~ —0.4 0.3 (i.e. z ~ 0.3 — 0.4), which can 

give some significant effects on ISW or WL. The growth index / obtained from general 7 (dotted 
line) is also consistent with the exact solution (solid one). The dashed line in / deviating from 
the exact solution is the one obtained when one use the constant 7 = 70 = 0.6. The only problem 
in this model is that it violates the solar system limit by F^/\/F\q ~ —0.7. If one release this 
limit which might be plausible in cosmological scale, then we may have a very interesting STG 
model which mimics exactly A model for the background evolution but shows the totally different 
behavior for the growth history. When we consider 70 < 0.56, the value of 4>' becomes imaginary 
at the present. Thus, STG with simple form of Z(cf>) = 1 may not be consistent with small value 
of 70 (i.e. negative 7 ). The shape of F(n)/Fo is determined by the sign of 7„. If 7„ is positive 
(negative), then F(n)/Fo shows the convex (concave) shape with the minimum (maximum) as 1. 
Thus, the value of F(n) is bigger than Fq in the past for the viable STG. 

Second, we probe the models which show the similar background evolutions as the some 
coupled quintessence models (cJo,u} a ) = (—1.0, 1.0). In these models, oj = during the tracking 
region and it approaches to —1 at present [83, 84]. Again we find 70 ^ 0.56 for j a = to 
get F(0) = 1. However, these values of STG model can not pass the solar system test because 
(F :4 ,/y/F)\ = -0.17. STG with other set of parameters (70,7a) ^ (0.562,0.006) gives F(0) = 1.00 
and (Fa/\/F~)\q = —0.02. Also <p' well behaves in this case. Thus, this model satisfy all known 
observational constraints. We show the evolutions of physical quantities of this in Fig. 2. Again oj 
obtained from Eq. (4.6) is exactly same as CPL oj with ojq = —1.0 and oj a = 1.0. F(n) decreases 
during the evolution which is opposite to the previous case even though j a > 0. However, this 
is only possible when 7 ~ 0. The growth index / obtained from general 7 (dotted line) is 
also consistent with the exact solution (solid line). Again / with constant value of 70 = 0.562 is 
deviated from the exact one. For (70,7a) — (0.6, 0.126), F(u)/Fq increases to 1.1 around n ~ —0.4 
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-2.0 -1.5 -1.0 -0.5 0.0 -2.0 -1.5 -1.0 -0.5 0.0 -2.0 -1.5 -1.0 -0.5 0.0 

n n n 

Figure 1: a) In the first row, we show the evolutions of u, F/Fq, and F s/yF (from left to right) 
when (—1.0,0) and (0.6,0.126) for (ojQ,oj a ) and (70, 7a), respectively, b) The evolutions of 4>, U, 
and / for the same values of u and 7 (from left to right) in the second row. 

(i.e. z ~ 0.4). This case also violates the solar system limit by F^/^/FIq ~ —0.45. 

In Table 1, we show the behaviors of physical quantities for various models with some parameter 
sets. There have been various quintessence models [85] but we just investigate the inverse power- 
law potentials (</>") [86, 87] with a = 1 which is well approximated by ojq = —0.74 and oj a = 0.07. 
STG for the same values of ujq and uj a with 70 = 0.57 is hardly distinguished from the 
quintessence model because / of both models are quite similar to each other. F(n)/FQ reaches 
to the maximum value 1.015 and then approaches to 1 at early time. No known observation can 
distinguish this from cfr^ 1 quintessence model. If (70,7a) are (0.6,0.08) for the same lj, then the 
maximum value of F(n)/Fo is 1.15 which can give some effects on ISW and WL. Also the time 
varying growth index parameter can be distinguished from the constant one. We also investigate 
the model inspired by the supergravity (SUGRA) [88]. The background evolution of SUGRA can 
be parameterized by (u>o, uj a ) = (—0.92, —0.08). Again, only STG with 70 = 0.6 with 7 a = 0.11 can 
give some significant deviations from SUGRA for the growth history. We also check the phantom 
crossing models in this table. Phantom crossing I means the model with uj > — 1 in the past 
becomes uj < — 1 at present. Phantom crossing II indicates the other model where oj < — 1 in the 
past and oo > —1 at present. The results do not change much for these cases. From this table, we 
find that only STG models with the high values of 70 can be distinguished from the other models. 
However, all of high values of 70 STG models violate the solar system test. 
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-2.0 -IS -1.0 -0.5 0.0 -2.0 -1.5 -1.0 -0.5 0.0 -2.0 -1.5 -1.0 -0.5 0.0 

n n n 

Figure 2: a) In the first row, we show the evolutions of u, F/Fq, and F^/\/F (from left to right) 
when (—1.0,1.0) and (0.562,0.006) for (wo,w a ) and (70,70); respectively, b) The evolutions of 4>, 
U, and / for the same values of u> and 7 (from left to right) in the second row. 



Models 






7o 


7a 


F(n)/F 


F,4>/Vf\ 


0' 


V(<t>) oc <j)- 1 


-0.74 


0.07 


0.56 


-0.018 


U min 0.986 


0.018 


fine 


0.57 





f| max 1.015 


-0.09 


fine 


0.6 


0.08 


f| max 1.150 


-0.48 


fine 


SUGRA 


-0.92 


-0.08 


0.56 


-0.016 


max 1.01 


0.00 


z < 1.3 


0.563 





H max 1.02 


-0.18 


fine 


0.6 


0.11 


H max 1.14 


-0.70 


fine 


Phantom 
Crossing 
I 


-1.1 


0.3 


0.53 


-0.09 


U min 0.93 


none 


imaginary 


0.557 





p| max 1.01 


-0.24 


fine 


0.6 


0.143 


p| max 1.15 


-0.68 


fine 


Phantom 
Crossing 
II 


-0.8 


-0.3 


0.55 


-0.049 


U min 0.97 


none 


imaginary 


0.568 





p| max 1.03 


-0.18 


z < 3 


0.6 


0.09 


p| max 1.13 


-0.69 


z < 7 



Table 1: Physical quantities for various models. U (P|) means the concave (convex) shape of 
F(ji)/Fq. Also /* indicates the monotonic increase of F/Fq as n increases. 4>' becomes imaginary 
after the maximum z for certain parameter set. 
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-2.0 -1.5 -1.0 -0.5 0.0 -2.0 -1.5 -1.0 -0.5 0.0 -2.0 -1.5 -1.0 -0.5 0.0 

n n n 

Figure 3: a) In the first row, we show the evolutions of to, F/Fq, and F s/vF (from left to right) 
for STG model with (70,7a) = (11/16,0.32). b) The evolutions of <fi, U, and / for the same values 
of to and 7 (from left to right) in the second row. 

4.2 Comparison with DGP model 

DGP model [30, 31, 32] is one of the well studied modified gravity models. The background 
evolution of this model can be parameterized by (tOQ,to a ) = (—0.78,-0.32) with 70 = 11/16 [17]. 
In Ref. [21], 70 = 4/7 at high z. We show the evolutions of physical quantities of STG with 
(7o>7a) = (11/16,0.32) in Fig. 3. Again except the solar system constraints, this model shows 
the interesting and different behaviors compared to DGP model = 0) for the growth history. 
F !( f,/F\o becomes very big negative value and this is due to the fact that (j)' approaches to at 
present which is used in F $ = F'/(j>'. Usually, / is smaller than 1 in general DE and MG models. 
However, we can produce / > 1 for 70 < 0.5 with non-vanishing j a models even though cf>' is 
imaginary. If observations do confirm this fact (/ > 1), then this kind STG might be a useful 
candidates to explain it. 

We show the current observational limits on / and 7 in Table 2. 7 value is obtained from / obs 
when we assume A model. The current observations may be consistent with A model except the 
Lyman-a forest result [95]. However, there exist huge errors in observations and it might be too 
early to extract any physical properties for specific models. 
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f2 m o 


jobs 


y>bs 


Ref 


0.15 


-0.14 


0.3 


0.51 ±0.11 


0-72™ 


2dFGRS [89, 90] 


0.32 


-0.28 


0.26 
±0.02 


0.6541S 


U -" JZ — 0.31 


SDSS R= 10 - 50/t- 1 Mpc [91] 


0.641±»;J«i 


o ^+ u - 2y 

u - oo -0.32 


R= 2 - bOh' 1 Mpc 


0.35 


-0.3 


0.3 


0.70 ± 0.18 


0.54+}$ 


SDSS [92] 


0.55 


-0.44 


0.3 


0.75 ±0.18 


0.59^ 


2dF-SDSS [93] 


1.4 


-0.88 


0.3 


0.90 ± 0.24 


0.68+f^ 


2dF-SDSS [94] 


3.0 


-1.39 


0.3 


1.46 ±0.29 


-10.6^ 


Ly-a (SDSS) [95] 



Table 2: is the corresponding redshift for each observation. 7 is the derived quantities from 
the observational value of / obs when we assume A model (u = —1). 

5 Conclusions 

Scalar-tensor gravities theories can produce many possible background evolutions which mimic 
dark energy models and other modified gravity models. We need to consider the growth of the 
linear matter perturbation to distinguish between models. However, STG models are strongly 
limited by the solar system constraint when we normalize the kinetic energy term Z(4>) = 1. The 
main reason for this is that (/>' becomes singular with this normalization. Thus, we may need to 
investigate the STG models with general Z(4>) in order to distinguish STG model with others. 

When we allow the time variation of the growth index parameter 7, usually the negative 
value of 7 a models have the singular problem of <j)' . Some models with the positive 7 a have the 
interesting features like large enough F/Fq values at early epoch while mimic the dark energy 
models background evolution. However, these cases violate the solar system test and again this 
might be able to be cured when we consider STG models. If j a is positive (negative), then F(ti)/Fq 
shows the convex (concave) shape with the minimum (maximum) as 1. Thus, the value of F(n) is 
bigger than Fq in the past for the viable STG. There can be exception for this case, when j a ~ 
but positive. 

The main conclusion is that the viable STG models with Z{4>) = 1 are not distinguishable 
from dark energy models or other modified gravity models when we strongly limit the solar system 
constraint. 

A Appendix 

Although the functional form of the growth factor given in Ref. [17] is simple and useful, its usage 
should be limited for certain values of u. If ojq and (or) u> a are (is) big, the form loses the accuracy. 
We show this in Fig. 4. We show the errors in 5 for the different values of ojq = —1.0 (solid), —0.8 
(dashed), and —0.6 (dotted) as a function of ui a in the left panel of Fig. 4. Except for oj$ = —1.0, 
the errors increase rapidly as u a increases. We also probe the errors for the different values of 
u a = 0.2 (solid), 0.4 (dashed), and 0.6 (dotted) as a function of ujq in the right panel of Fig. 4. 
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O B— - ■ ■ ■ 

-0.4 -0.2 0.0 0.2 0.4 -1.2 -1.1 -1.0 -0.9 -0.8 -0.7 -0.6 

<jj a CJq 

Figure 4: a) Errors in the growth factor <5 obtained from Eq. (1.1) compared to the exact 
(numerical) solution as a function of oj a for the different values of ojq = —1.0, —0.8, and —0.6 
(from bottom to top), b) Errors as a function of ojq for the different values of u a = 0.2, 0.4, and 
0.6 (from left to right). 

The errors are more than 1 % in many cases. Thus, one should not rely on the functional form 
for general case. 



Also this functional form does not depend on f2 m o. Thus, if one study some model with small 
or large enough values of Sl m o, then one should not rely on this functional form. Instead one can 
improve 7 by using the exact solution of 5 as shown in Refs. [14, 15, 16] 
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in[n 



mOj 



(A.l) 



where Q 



±jgf and w cff = u + f 



The initial values of <j)' and U{4>) are determined if we know the exact values of (wo,w a ), 
(7o, 7a), and f2^, which might be possible up to certain level in future observations. We use Eqs. 
(2.7), (4.6), (4.2), (4.8), and (4.9) to obtain ^ and U((j> ). 



2F 
Up 
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3F 



Hq ' 3Fq 



- ^mo) 
^mo) , 



(A.2) 
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